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Abstract 

Using an isospin- and momentum-dependent transport model we study the emission of free nu- 
cleons and the nuclear radial flow in central heavy-ion collisions induced by high energy radioactive 
beams. The midrapidity neutron/proton ratio and its transverse momentum dependence are found 
very sensitive to the high density behavior of nuclear symmetry energy. The nuclear radial flow, 
however, depends only weakly on the symmetry energy. 
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The density dependence of nuclear symmetry energy E 
derstanding many interesting questions in astrophysics [l[ |_ 

of radioactive nuclei 0, Q, si, Q] and many observables of heavy- ion reactions 0, rJ 11, fl3 | . 




is^ery important for un- 
the novel structures 



However, it is still very poorly known both theoretically and experimentally, especially at 
supranormal densities, see e.g., refs.P, 13]. Fortunately, nuclear reactions induced by ra- 
dioactive beams provide a unique opportunity to pin down the symmetry energy in a broad 
density range Q]. Moreover, they also allow us to explore the role of isospin degree of free- 
dom in various phenomena and the dynamics of nuclear reactions. In particular, central 
heavy-ion reactions induced by high energy radioactive beams up to 400 MeV/nucleon at 
the planned Rare Isotope Accelerator (RIA) will enable us to probe the high density behav- 
ior of nuclear symmetry energy. While a number of experimental observables constraining 
the E sym (p) at subnormal densities have been identified recently, see e.g., refs. P. ll(l \v\ for 
reviews, very few observables sensitive to the high density behavior of E sym (p), except the 
7y~/ti + ratio jlfii llfi|. are currently known. In this work, we examine midrapidity nucleons 
and the radial flow in central heavy-ion reactions induced by high energy radioactive beams 
as potential probes of the symmetry energy at high densities. 

Our study is based on the isospin- and momentum-dependent transport model 
IBUU04 117| . In this model a parameter x was introduced in the single nucleon poten- 
tial to mimic various density dependences of the symmetry energy as predicted by different 
microscopic many-body theories^]. Shown in Fig. 1 are examples of the symmetry energy 
with the parameter x going from 1 (softer) to -2 (stiffer). The recent data from NSCL/MSU 
on isospin diffusion in mid-central 124 Sn + 112 Sn reactions at 50 MeV/nucleon was used to 
constrain the symmetry energy at subnormal densities Within the IBUU04 model, a 
symmetry energy E sym (p) 31.6(p// ?n) 1 ,05 corresponding to x = — 1 was found to reproduce 



the isospin diffusion data very nicely 20(. In the present study we use the same model and 
parameter set to study nucleon emissions and the radial flow for central 132 Sn + 124 Sn reac- 
tions at a beam energy of 400 MeV/nucleon. In this reaction a compressed matter of about 
twice normal nuclear matter density is formed and lasts for about 15 fm/c[16]. This par- 
ticular reaction system will be studied experimentally using the Time Projection Chamber 
with beams from the fast fragmentation line at RIA. 

Shown in Fig. 2 are the rapidity distributions of free nucleons identified as those having 
local baryon densities less than po/8 at 30 fm/c when the dynamical freeze-out for particle 
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emission is reached. It is seen that mid-rapidity nucleons are very sensitive to the symmetry 
energy. With the symmetry energy changing from being softer (x=l) to stiffer (x=-2) 
more (less) neutrons (protons) are emitted at midrapidity. This is a direct consequence of 
the higher repulsive (attractive) symmetry potential for neutrons (protons) with the stiffer 
symmetry energy at supranormal densities. Moreover, neutrons appear to be more sensitive 
to the variation of the symmetry energy than protons. This is because for protons the 
repulsive Coulomb potential works against the symmetry potential. It is also interesting 
to point out that the effects of the symmetry potential on nucleon emissions observed here 
are just opposite to those at intermediate energies around 50 MeV/nucleon where a softer 
symmetry energy causes more (less) emissions of neutrons (protons) 21]. This is because of 
the crossover of the symmetry energy functions with different density dependences at p as 
shown in Fig. 1. Thus it is very important to carry out excitation function studies of particle 
emissions in order to map out the symmetry energy in a broad range of density. 

The degree of isospin equilibrium or translucency can be measured by the rapidity dis- 
tribution of nucleon isospin asymmetry £yw = (N n — N p )/(N n + N p ) where N n (N p ) is 
the multiplicity of free neutrons (protons) [22] . Although it might be difficult to measure 
directly St ree because it requires the detection of neutrons, similar information can be ex- 
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tracted from ratios o 
coalescence model 123, 



Jight clusters, such as, 3 H/ 3 He, as demonstrated recently within a 



24j . Shown in Fig. 3 are the rapidity distributions of Sf ree with (upper 



window) and without (lower window) the Coulomb potential. It is interesting to see that 
the Sf ree at midrapidity is particularly sensitive to the symmetry energy. As the parameter 
x increases from —2 to 1 the <5/ ree at midrapidity decreases by about a factor of 2. Moreover, 
the forward-backward asymmetric rapidity distributions of £/ ree with all four x parameters 
indicates the apparent nuclear translucency during the reaction. Comparing the results 
shown in the upper and lower windows, one sees that the Coulomb potential is to reduce 
the value of Sf ree , while the sensitivity to the symmetry energy remains about the same. 

Concentrating on nucleons at midrapidity, we examine in Fig. 4 the differential neu- 
tron/proton ratio dN n /dNp as a function of transverse momentum p t . The solid line in 
the figure is the average (n/p) sys ratio of the reaction system. In the low (high) p t part 
the dN n /dN p is significantly higher (lower) than the (n/p) sys . Moreover, the low p t part of 
the dN n /dN p ratio is more sensitive to the symmetry energy than the high p t part. Both 
observations are due to the momentum dependence of the symmetry potential. As shown 
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earlier by one of us the symmetry potential decreases with increasing nucleon momentum in 
agreement with the Lane potential extracted from nucleon-nucleus scattering data^l l^ . 
Of course, the symmetry potential is also density dependent. Generally, high p t particles 
are more likely to come from high density regions where the symmetry potential is stronger. 
However, if effects of the momentum dependence are stronger than those of the density de- 
pendence, the low p t particles are then more sensitive to the symmetry potential as observed 
here. 

It is also seen from Fig. 4 that the dN n / dN p ratios with different x parameters are almost 
parallel with each other except in the high p t region. It is well known within the blast wave 
picture of central havey-ion reactions that the inverse slope of the particle spectra, i.e., the 
apparent temperature, reflects the combined effects of the true temperature and the radial 
flow velocity of the compressed matter in the exploding fireball[26]. Our observation here 
about the slopes of the dN n /dN p ratios indicates that the symmetry energy has a weak 
influence on the overall thermalization and the radial flow. To be more specific we study 
in the following the isospin dependence of radial flow. It is worth noting that to our best 
knowledge this kind of study has never been carried out before. In anticipation of the coming 
experiments at RIA, it is important to identify what role the isospin degree of freedom may 
play in the radial flow which is in its own right an interesting phenomenon. We extract the 
average radial flow velocity at radius r from 

1 - 

where N(r) is the number of particles including bound ones in the spherical shell between 
radius r and r + dr, pi, r$ and are the momentum, coordinate and energy of the particle 
i. As an example, we show in Fig 5 the radial flow velocities of neutrons and protons and 
their densities as functions of radius r measured from the center of mass of the reaction with 
x = — 1. It is seen that the radial flow velocity increases with the radius r similar to the 
Hubble expansion of the universe. This observation is consistent with other studies on the 
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281 ] . Protons, especially at larger radii, 



radial flow in heavy- ion reactions, see, e.g., refs. 
are flowing with a slightly higher velocity due to the Coulomb repulsion, although they have 
about the same density profile as neutrons. To explore effects of the symmetry energy on 
the radial flow we show in Fig. 6 the position-averaged nucleon radial flow velocity < (3 > 
with and without the Coulomb potential as a function of the parameter x. The value of 



< (3 > can be readily extracted experimentally from the particle spectra using the Siemens- 
Rasmussen formula 2fJ. In our calculations it is obtained from eq^by setting r = and 
dr = oo. Formally, < (3 > is related to (3{r) via 



where N tota i is the total number of neutrons or protons. Without the Coulomb potential, 
the radial flow of neutrons is faster than that of protons because of the repulsive (attractive) 
symmetry potential for neutrons (protons). The difference in < (3 > for neutrons and 
protons is the largest for the stiff est symmetry energy considered here as one expects. As the 
symmetry energy becomes softer the difference disappears gradually. However, the overall 
effect of the symmetry energy on the radial flow is small, even for the stiffest symmetry 
energy with x = —2 the effect is only about 4%. This is because the pressure of the fireball 
is dominated by the kinetic contribution. Moreover, the compressional contribution to the 
pressure is overwhelmingly dominated by the isoscalar interactions. For protons, however, as 
the Coulomb potential is turned on, because of its dominance over the symmetry potential 
the whole picture changes. In fact, the Coulomb potential almost cancels out the effect of 
the symmetry potential at x — —2. As the symmetry energy becomes softer, the radial 
flow for protons becomes higher than that for neutrons. With the softest symmetry energy 
considered here, i.e, with the parameter x = 1, the Coulomb potential results in a higher 
radial flow velocity for protons than neutrons by about 3%. The radial flow for neutrons 
is only weakly affected through a secondary effect because the reaction dynamics is slightly 
modified by the Coulomb interaction. 

In summary, using an isospin- and momentum-dependent transport model we have stud- 
ied nucleon emissions and the nuclear radial flow as potential probes of the symmetry energy 
at high densities. We found that the isospin asymmetry of midrapidity nucleons and its 
transverse momentum dependence are very sensitive to the high density behavior of nuclear 
symmetry energy. These observables together with the tt~ /tt + ratio reported earlier will by 
very useful for studying the equation of state of dense neutron-rich matter at RIA. We have 
also studied the isospin dependence of nuclear radial flow. Our findings clearly indicate that 
the isospin degree of freedom plays an important role in the nuclear radial flow through the 
competition between the Coulomb and the symmetry potential. However, the net effect of 
the symmetry energy on the radial flow is small. 
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FIG. 1: (color online) Nuclear symmetry energy as a function of density with different x parameter. 
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beam'c.m. 

FIG. 2: (color online) Rapidity distributions of free neutrons and protons in the reaction of 
132 Sn + 124 Sn at Eb eam /A = 400 MeV and an impact parameter of 1 fm. 
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FIG. 3: (color online) Rapidity distribution of isospin asymmetry of free nucleons in the same 
reaction as in Fig. 2 with (upper window) and without (lower window) the Coulomb potential. 
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FIG. 4: (color online) Transverse momentum dependence of the differential neutron/proton ratio 
at mid-rapidity in the same reaction as in Fig. 2. 
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FIG. 5: (color online) Nuclear radial flow velocity and density as a function of position for the 
same reaction as in Fig. 2. 
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FIG. 6: (color online) Symmetry energy dependence of the position-averaged radial flow velocity 
for the same reaction as in Fig. 2. 
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